Orbital ordering in charge transfer insulators 
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We discuss a new mechanism of orbital ordering, which in charge transfer insulators is more 
important than the usual exchange interactions and which can make the very type of the ground 
state of a charge transfer insulator, i.e. its orbital and magnetic ordering, different from that of 
a Mott-Hubbard insulator. This purely electronic mechanism allows us to explain why orbitals in 
Jahn- Teller materials typically order at higher temperatures than spins, and to understand the type 
of orbital ordering in a number of materials, e.g. K2CUF4, without invoking the electron-lattice 
interaction. 
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Since the seminal Zaanen, Sawatzky, and Allen paper 
0, it is widely accepted that strongly correlated insula- 
tors can be divided into two classes according to the na- 
ture of their gap. The Mott-Hubbard insulators (MHI) 
are characterized by the lowest-energy excitations of the 
type d n +d n — > d n_1 +d n+1 , obtained by transferring a d- 
electron from one transition metal (TM) ion to another, 
whereas in charge transfer insulators (CTI) the corre- 
sponding excitations are holes on ligand sites (oxygen 
ions in the case of TM oxides): d n +p 6 -> d n+1 +p 5 . MHI 
are usually described by an effective Hubbard model, in 
which only the d-electron states are considered explicitly. 
On the other hand, CTI are described by models that in- 
clude also the relevant oxygen p-states, such as the dp (or 
three-band Hubbard) model of high-T c superconductors 
Q. The value of the gap is determined by the smallest 
of the two parameters: the energy of electron transfer 
from ligand to TM sites A and the Hubbard repulsion 
on TM sites Ud- For CTI (usually, the oxides of heavier 
3d metals) A < U d , while for MHI U d < A. 

The distinction between these two types of insulators 
becomes apparent upon doping and may play an impor- 
tant role in high-T c superconductors |2, L|- though ac- 
cording to other point of view, there is no qualitative 
difference between doped CTI and MHI |^|. In the ab- 
sence of doping, both CTI and MHI with one electron or 
hole per TM site are S = \ Heisenberg antiferromagnets. 

The situation turns out to be quite different in Jahn- 
Teller systems, in which d-electrons, in addition to spins, 
have also orbital degrees of freedom. In this Letter we 
show that in insulators with orbital degeneracy not only 
the excitation spectrum, but also an orbital and magnetic 
ordering in the ground state may strongly depend on the 
ratio A/Ud- The difference between the ground states of 
CTI and MHI results from a peculiar exchange interac- 
tion between orbitals on nearest-neighbor d-sites, which 
may be called 'the orbital Casimir force', as it originates 
from the vacuum energy of charge fluctuations due to 
the hopping of electrons between the ligand and TM ions 
(whereas the usual exchange involves the electron trans- 



fer from one TM site to another). We recently noted 
the importance of such interaction in geometrically frus- 
trated systems with the 90° metal-oxygen-metal bonds, 
e.g. NaNi02, where it determines the type of orbital 
ordering p| ■ Here we show that this orbital exchange of 
electronic origin is generally present in all TM oxides and, 
in particular, in materials with the 180° metal-oxygen- 
metal bonds, like perovskites ABO3 and their layered 
analogs A2BO4, and it is not included in the widely used 
exchange model jy, [j] (sometimes called Kugel-Khomskii 
(KK) model). That model is derived from an effective 
Hubbard model describing only the degenerate d-electron 
orbitals and which, strictly speaking, only applies to MHI 
with A -C Ud- The new orbital interaction obtained be- 
low helps to resolve a difficulty of the KK model, in which 
an orbital ordering is strongly coupled to a magnetic one, 
implyingclose temperatures of orbital and magnetic or- 
derings 7], whereas in reality orbitals usually order at 
higher temperatures than spins. 

Exchange interactions: We first calculate an exchange 
Hamiltonian for two TM ions with one hole on the doubly 
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The two orthogo- 
3z 2 — r 2 ) and \x 2 — y 2 ) may be described 
by, respectively, the up and down states of the isospin- 
\ operator T z 0, 0- It is convenient to introduce an 
overcomplete set of operators d iaa on each TM site i an- 
nihilating a hole in the state |3a 2 — r 2 ) (a = x, y, z) with 
the spin projection a. These states are hybridized with 
the p a orbitals on two neighboring oxygen sites i ± ^ 
(the corresponding oxygen hole operators are denoted by 
Pi+SL iCT ). The hole states on TM sites are described by the 
standard atomic Hamiltonian, which includes the on-site 
repulsion Ud and the Hund's rule exchange coupling Jd 
for holes on two orthogonal orbitals. The Hamiltonian of 
oxygen holes hybridized with the d-states is 
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FIG. 1: Diagrams describing exchange interactions mediated 
by two virtual oxygen holes. 



where the first term is the kinetic energy of holes in the 
oxygen p-band, the second term is the Hubbard repulsion 
on oxygen sites, and the last term describes the hopping 
between oxygen and TM sites. 

Integrating the oxygen degrees of freedom out, one ob- 
tains an effective interaction between (i-electrons in the 
form of an expansion in powers of (t p d/A) 2 . The lowest 
order term provides a hopping between TM sites, which 
together with the Coulomb and Hund's rule interactions 
between (i-electrons gives an effective dd-model. For dis- 
persionless oxygen holes, H k i n = & Hi. a ,<jP\+^<jPi+ ! i<n 
the effective dd-hopping occurs between nearest-neighbor 
TM sites, and its amplitude is tdd — i^ d / A. 

The next term in this expansion (oc t p d ) gives an or- 
bital and spin exchange mediated by two oxygen holes. 
The corresponding diagrams are shown in Fig^ where in 
the diagrams (a) and (b) two holes do not interact with 
each other in the virtual state, while the black squares 
on diagrams (c) and (d) correspond to one or more Hub- 
bard interactions on oxygen sites. Each diagram gives 
rise to a particular form of spin and orbital exchange. 
Thus, if the holes on the TM sites i and j exchange 
places [Figda)], the spin exchange is described by oper- 
ator P ex — +1 if the total spin of the two holes S = 1, 
and —1 for 5 = 0. This operator can also be written 
in the form P ex = 1 - 2P , where P = 1/4 - (SjSj) is 
the projector on the S — state. If, on the other hand, 
the holes return back to their starting points, the cor- 
responding 'spin factor' is —1. The interacting holes in 
the diagrams (c) and (d) must have opposite spin pro- 
jections, so the spin exchange in that case is described 
by S^Sj +S^~ Sj' [see Fig^c)], while if the spin projec- 
tions remain unchanged [see Fig^d)], the spin operator 
is —2 (1/4 — SfSj). Since the diagrams (c) and (d) only 
differ by the spin factors, their sum is proportional to 
— 2Pq, corresponding to the fact that the on-site Hub- 
bard interaction between two holes is only nonzero when 
their spins add to the total spin 0. 

While the spin-dependence of the exchange interac- 
tions, discussed above, is very general, the orbital struc- 



ture depends on the dispersion of holes and on relative 
positions of the sites i and j. In particular, for disper- 
sionless holes the exchange in order t* d only occurs be- 
tween two nearest neighbor TM sites, e.g. i and i + a, 
and the diagrams (a), (c), and (d) are proportional to 



T-ict — ( 2 + £f) (2 
1 V3 



T x ry 



where 



-T z - 



V3 



rpX JZ 



(2) 



and the operator (| + if) = d\ a di a is the projector 
on the state |3a 2 — r 2 ) on the site i. The contri- 
bution of the diagram in Fig^b) is proportional to 
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Combining terms containing Pq with the exchange in- 
teractions appearing in the dd- model in the order tp d , we 
obtain 
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where the last term describes the Hund's rule coupling on 
TM sites and P x = 1 -P = 3/4+ (SjSj) is the projector 
on the S = 1 state. In addition to the coupled spin and 
orbital interactions, there exists also a spin-independent 
orbital exchange 
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The first term is the spin-independent contribution of 
the processes shown in Figs^a) and (b), while the sec- 
ond one results from the dd-model. It is important to 
note that in the order t^ d it is impossible to exchange 
two holes in such a way that they never occupy the same 
oxygen site (the corresponding terms in the diagrams (a) 
and (b) cancel each other). The orbital exchange inter- 
action containing only powers of A is actually due to the 
virtual hopping of holes to all neighboring oxygen sites 
and back [Fig^b)], except for the oxygen bridging the 
two neighboring TM sites. In the limit A > Ud Eqs. © 
and reduce to the KK model. 

The term in Ht containing only A in denominator is 
similar to the pure orbital exchange for the 90° metal- 
oxygen- metal bonds 0, and its origin can be better un- 
derstood by noting that it is the only term that does not 
vanish in the limit of infinite on-site Hubbard repulsion 
Up and Ud- It represents the energy of the derealization 
of holes over the octahedra of the ligand ions, surround- 
ing the TM sites, or more correctly, the increase of that 
energy due to the fact that the neighboring octahedra 
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FIG. 2: Favorable [(a) and (b)] and unfavorable (c) occupa- 
tions of hole orbitals of neighboring TM ions in CTI, in which 
the exchange is dominated by Eq.(0J (the terms oc Uf 1 sta- 
bilize the state (a)). 



share one ligand ion and for infinite U p the intermediate 
two-hole states on the shared ion are forbidden. This 
energy increase depends on orbital states of the neigh- 
boring holes, which gives rise to the orbital exchange 
with the positive sign. Thus the origin of this orbital 
interaction is similar to that of the Casimir effect Q , re- 
sulting from the dependence of the energy of zero-point 
fluctuations of electromagnetic field on positions of two 
metallic plates. It is also important that this orbital ex- 
change has positive sign, while all other terms in Eqs. © 
and (0} are negative (since they describe the second-order 
corrections to the ground state energy in the effective (id- 
hoping). The result of this new orbital interaction can 
be formulated as the following "rule of thumb" (applied 
also to materials with one electron per TM site) : to avoid 
the energetically unfavorable occupation of the bridging 
oxygen by two holes, the d-holes should occupy orbitals 
directed away from the common oxygen either on both 
neighboring TM sites [see Figl^b)] or, at least, on one 
site [see FigE^a)]. 

The importance of the orbital exchange Eq.|0J) for an 
orbital ordering may be understood already by consider- 
ing the minimal-energy states of a pair of holes on the 
neighboring TM sites i and i + z. In MHI (A > U d ) the 
exchange is dominated by Hst © and the last term in 
Eq.Q. The AFM exchange between spins of the holes, 
described by this term, is the strongest for I z = i on 
both sites, i.e., the lowest energy is achieved if both holes 
occupy the |3z 2 — r 2 ) state, as shown in FigGfc). The 
ground state of the pair is a spin singlet. On the other 
hand, in CTI (A < Ud) the strongest interaction is the 
orbital exchange Ht (see Eq.Q), and the state shown in 
FiglUc) has the highest orbital energy, since both hole or- 
bitals are oriented towards the common ligand ion. The 
more energetically favorable states in CTI will be those, 
in which at least one hole is in the \x 2 — y 2 ) state (see 
Fig|2K and b). Furthermore, both the last term in Eq.QJ 
and the relatively weak FM spin exchange due to the 



Hund's rule coupling select the state shown in Fig^ with 
the total spin of the pair S = 1. Thus, the ground state 
of the pair in the CT regime is different from that in the 
MH regime. 

Consider now two neighboring TM sites with one e g 
electron, as e.g. Mn 3+ (i^ej), Cr 2+ (^l ff e^), and the 
low-spin Ni 3+ (t^eg)). In the KK model the systems 
with 1 electron/site and 1 hole/site order in exactly the 
same way, since the Hubbard model, from which the KK 
exchange Hamiltonian is derived, has the electron-hole 
symmetry 0; consequently, the term Hst © has the 
same form both for electrons and holes. On the other 
hand, in the dp-model this symmetry is lost: the new 
orbital exchange interaction for one-electron case has the 
form [cf. Eq.Ql] 




In this case for both MHI and CTI the energy of a pair 
of TM electrons on the sites i and i + z is minimal if both 
electrons occupy the d 3z 2_ r 2 orbital, Figl^c), and their 
spins form the singlet state. 

Note that the electron-phonon (Jahn- Teller) interac- 
tion also gives rise to a purely orbital exchange of the 
form Hjt oc Ylia^i^i+a^ wmcn f° r both the electron 
and hole cases favors the antifcrro-orbital (AFO) spin- 
triplet state shown in FigEJa). Thus the JT interaction 
favors the same state as the novel orbital exchange for 
one-hole systems, but a different one for one-electron sys- 
tems. This, in principle, gives a possibility to experimen- 
tally discriminate between these two mechanisms and to 
study their relative importance. The Cu 2+ pairs (one 
hole) in K2Z11F4 were found to have the triplet ground 
state |9j, showing that the orbital exchange Eq.(0} and/or 
JT interaction are in this case stronger than the KK ex- 
change. We are not aware of similar data for pairs of 
one-electron ions. 

Orbital and magnetic ordering: Equations allow 
us to consider types of ordering in different situations. 
Thus, e.g., the ground state of a chain of TM ions in 
the z-direction with 1 hole/site is AFO with alternating 
d^ z 2_ r 2 and d x 2_ y 2 orbitals and FM for CTI, as that is 
favored by the sum of Eqs. (J3J) and while for MHI, 
in which Hst dominates, we obtain the ferro-orbital or- 
dering with <i 3z 2_ r 2 orbitals occupied on all sites, and 
the AFM spin exchange 01- The latter state also mini- 
mizes the exchange energy of one-dimensional CTI with 
1 electron/site. 

Similarly e.g. for CTI (Ud — > 00) with a square lattice 
and 1 hole/site, the orbital exchange interaction Eq.QJ 
has the form 

Hr=+^j:m +x +iii! +y -i!). (6) 
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The term bi-quadratic in isospin operators favors an an- 
tiparallel orientation of neighboring isospins, while the 
linear term in Eq.© acts as 'the orbital field' in the z 
direction, favoring the canted orbital state with the cant- 
ing angle ±60°, which corresponds to the two-sublatticc 
orbital ordering with the d z 2_ x 2 and d z 2_ y 2 hole states. 
Since these two neighboring orbitals do not overlap, the 
Hund's rule coupling (the last term in Eq.(PJ) results in 
the FM exchange between the spins of the holes. This is 
precisely the orbital and magnetic structure of the lay- 
ered compound K2CUCI4. For MHI, however, we would 
get from Eq.@ a different ground state: ferro-orbitally 
ordered with I z — — \ (d^.2 _,.2 -orbitals) and AFM spin 
ordering (see, however, Ref. Il4|) . Previously, the KK 
model was successfully used to predict the orbital struc- 
ture of that compound |lCj . assuming that spins order 
ferromagnetically; but as we mentioned above, such a 
state is not the ground state of the KK model. 

The case of a cubic lattice of TM ions is somewhat 
special for several reasons. First, using Yla^i ~ 0> we 
find that the exchange interactions have the electron-hole 
symmetry. In particular, both Eq.|0J and Eq.JSJl can be 
written in the form 




Second, the orbital interaction Eq.Q is frustrated 
[T3 | : apart from the AFO mean-field ground states with 
the average isospin vector (Tj) parallel (antiparallel) to 
an arbitrary vector in the (T x , T z )-plane, there are in- 
finitely many disordered ground states, obtained from the 
AFO states by performing the transformation I a — ► I a , 
I 13 JTj where a, (3, and 7 are three different spatial 
indexes, on an arbitrary set of planes transverse to the a 
axis. This ground state degeneracy is similar to the one 
of the orbital Hamiltonian of e g electrons on a triangular 
lattice [5j- The quantum orbital fluctuations select the 
AFO states I a = ±5, where a = x,y or z, which have 
the softest (two-dimensional) orbital excitations ^3] • F° r 
instance, for a = z, such a state corresponds to the oc- 
cupation of d 3z 2_ r 2 and d x 2_ y 2 orbitals on A and B sub- 
lattices. The corresponding magnetic ordering would be 
AFM in the xy planes and FM in the z direction. 

If we now would increase the ratio A/Ud, moving from 
CTI to MHI, the coupling between orbitals and spins 
would ultimately stabilize the canted AFO state with an 
AFM coupling along the chains in the z direction and a 
FM exchange between spins in the xy layers. The de- 
tailed treatment of this case, relevant e.g. for the per- 
ovskites KC11F3 and LaMn03, will be given elsewhere. 

In conclusion, we showed that the description and 
properties of orbital ordering in charge transfer insula- 
tors are much different from those in Mott-Hubbard in- 



sulators. In particular, even in perovskite CTI with the 
180°-exchange the orbital interactions are stronger than 
the spin exchange, so that an orbital ordering, in gen- 
eral, occurs at higher temperatures. In this respect the 
exchange mechanism of an orbital ordering in CTI resem- 
bles the ordering due to the JT (electron-phonon) inter- 
action, although is not identical to it. Another difference 
between the CTI and MHI is that the electron-hole sym- 
metry existing in the MH case is lost in the CT regime, so 
that the type of an orbital and spin ordering may be dif- 
ferent for one-electron and one-hole materials. The effec- 
tive exchange Hamiltonian for CTI allows us to describe 
the ordering in a number of such materials, without in- 
voking any extra factors. Most importantly, in orbitally 
degenerate insulators not only the excitation spectrum, 
but also the very type of the ground state may be dif- 
ferent in Mott-Hubbard and charge transfer regimes, in 
contrast to the non-degenerate case, in which it is the 
same. 
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